Mean field study of structural changes in Pt isotopes with the Gogny interaction. 
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The evolution of the nuclear shapes along the triaxial landscape is studied in the Pt isotopic chain 
using the selfconsistent Hartree-Fock-Bogoliubov approximation based on the Gogny interaction. 
In addition to the parametrization DIS, the new incarnations DIN and DIM of this force are 
also included in our analysis to asses to which extent the predictions are independent of details 
of the effective interaction. The considered range of neutron numbers 88 < A'^ < 126 includes 
prolate, triaxial, oblate and spherical ground state shapes and serves for a detailed comparison 
of the predictions obtained with the new sets DIN and DIM against the ones provided by the 
standard parametrization Gogny-DlS in a region of the nuclear landscape for which experimental 
and theoretical fingerprints of shape transitions have been found. Structural evolution along the Pt 
chain is discussed in terms of the deformation dependence of single particle energies. 

PACS numbers: 21.60.Jz, 27.70.-(-q, 27.80. +w 
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I. INTRODUCTION. 

Quadrupole collectivity is one of the most relevant fea- 
tures of nuclear structure [H, In this context, the 
theoretical understanding of the evolution of the nuclear 
shapes, and the structural changes associated with it, 
represent an active research field [3|-|9[. From the exper- 
imental side, low-lying spectroscopy is one of the most 
powerful sources of information about nuclear shapes 
and/or shape transitions since one can establish signa- 
tures correlating the excitation energies with the defor- 
mation properties 

Nowhere, however, the evolution of nuclear shapes is 
more documented and challenging than around the pro- 
ton shell closure Z = 82. For example, the neutron 
deficient Lead isotopes with neutron number iV w 104 
display three 0+ states within 1 MeV excitation energy 
[I?]. The very rich and challenging variety of nuclear 
shapes also extends to the neighboring Hg and Po iso- 
topes [l^l- In particular, it has been demonstrated [ll] 
that in the case of the Lead isotopes the decreasing trend 
observed in the binding energy difference 62p{Z,N) — 
E{Z -2,N) - 2E{Z,N) + E{Z + 2,N) for decreasing 
mass number A, can already be described quantitatively 
by mean field models in terms of deformed ground states 
of Hg and Po nuclei while the inclusion of the quadrupole 
correlation energy [191] brings the calculations even closer 
to experiment. From the experimental point of view [lo| . 
the neutron deficient Mercury isotopes exhibit deformed 
ground states while the situation is more involved in the 
case of Po nuclei (see, for example 20:], and references 
therein). 

A considerable effort, has also been devoted to charac- 
terize Pt nuclei [la [l3l. [2l| - [29| . In this case, several defor- 
mation regimes have been suggested. Previous theoreti- 
cal investigations 30. - .35] have found triaxial and oblate 
ground state shapes for the heaviest Pt isotopes while 
for the light ones a prolate deformed regime is predicted. 



From the experimental point of view [Hill, Hi the en- 
ergy ratio E4+/E2+ is almost 2.5 for Pt nuclei with neu- 
tron number 110 < A'' < 118 already pointing to 7 soft 
shapes. The role played by the 7 degree of freedom in 
Pt isotopes has also been stressed by the comparison of 
experimental and theoretical results performed in Ref. 
37] which shows, that good agreement can be obtained 
if triaxiality is taken into account. Further down, a tran- 
sition to a vibrational regime is suggested for i68-i72p|- 
by both the experimental data and their theoretical in- 
terpretation |22| . 

The shape evolution provided by the mean field frame- 
work ^], based on the most recent incarnations of the 
Gogny interaction jsl], in the isotopic chain i66-204p|- 
considered in the present study as a representative sam- 
ple of nuclei close to the Z = 82 proton shell closure for 
which prolate, triaxial, oblate and spherical ground state 
shapes are found. 

Nuclear shapes around the proton magic number Z = 
82 have been studied using a wide variety of theoretical 
models. Low-lying minima in both Pb and Hg isotopes, 
have been predicted within the framework of the Struti- 
nsky method [39l - l4l| . From a mean field perspective, the 
coexistence between different nuclear shapes in Pt, Hg 
and Pb nuclei, has been considered within the relativis- 
tic mean field approximation [4^ . Deformed ground 
states were predicted in the case of Pb isotopes as well as 
superdeformed ground states in Hg isotopes at variance 
with experiment. In order to cure this problem, a new 
parametrization of the relativistic mean field Lagrangian, 
called NLSC, was introduced in Ref. 
a new set called NL3* has been proposed in Ref. 



More recently, 

providing an improved description of the ground state 
properties of many nuclei like Pb isotopes. Studies based 
on Skyrme-like models have been reported [iol-lsoj while 
the shape coexistence in i82-i92pj-| ^^^^ analyzed in Ref. 
[5l| using the Hartree-Fock-Bogoliubov (HFB) approach 
based on the parametrization DIS ^] of the Gogny in- 
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teraction 38[. From a beyond mean field perspective, 
symmetry projected configuration mixing, based on both 
Skyrme [53, [s^ and Gogny [s^ energy density function- 
al, lias been successfully employed in this region of the 
nuclear chart establishing a firm ground to support the 
experimental evidences for rotational bands in the neu- 
tron deficient Pb isotopes built, on coexisting low- lying 
0+ states. On the other hand, evidences for 7 vibra- 
tions and shape evolution have also been considered in 
the nuclei i84-i90jjg jg^^ where a five-dimension collec- 
tive Hamiltonian was built with the help of constrained 
Gogny-DlS HFB calculations. 

Just below the proton magic number Z = 82, the nu- 
clei with A — 170 — 200 are particularly interesting be- 
cause, small islands of oblate deformation might be fa- 
vored energetically. Transitions from prolate to oblate 
shapes, as the number of neutrons increases, have been 
predicted in this mass region, using collective models 
p^3i |. phenomenological Woods-Saxon or Nilsson poten- 
tials jH, [13, relativistic mean field [soj - fBlj as well 
as non-relativistic deformed Hartree-Fock (HF) [1, [g^I 
and HFB calculations [8, 63] . Signatures for a transition 
from prolate to oblate ground states, as the number of 
neutrons increases from N = 110 to = 122, have been 
found in Hf, W and Os isotopes. In particular, such a 
transition was found [Ij to happen at = 116 — 118. 
Subsequently, the evolution of the ground state shapes 
along the triaxial landscape of several isotopes of Yb, Hf, 
W, Os and Pt has been studied in Ref. [63] within the 
framework of the mean field approximation based on both 
Skyrme and Gogny interactions. This region is also an 
active research field within the Interacting Boson Model 
(IBM) [64-72]. 

Taking into account that around the neutron mid-shell 
N=104 examples of coexisting configurations have been 
found [10.] , it is very interesting to study the propagation 
of the nuclear shapes in the Pt isotopic chain and the 
way it can be correlated with the details in the underly- 
ing single-particle levels as functions of the deformation 
parameters. For such a task, the mean field approxi- 
mation appears as a first tool incorporating, important 
correlations within the concept of spontaneous symme- 
try breaking and allowing a description of the evolution 
of shell structure and deformation all over the nuclear 
chart (see, for example, d, HE IzH and references 
therein) . 

In the present work, our study will be performed within 
the selfconsistent HFB framework based on the Gogny 
interaction [s^. In addition to DIS [s^j, which is still 
the most standard and thoroughly tested parametriza- 
tion, we also considered the two most recent parame- 
ter sets of the Gogny interaction, i.e., DIN [TJ] and 
DIM (75|]. To the best of our knowledge, the results 
to be discussed later on in this paper, are the first sys- 
tematic mean field study reported in the literature, us- 
ing both DIN and DIM parametrizations, to describe 
(mean field) ground state properties of Pt nuclei. The 
selected isotopes i66-204p^ cover almost the whole shell 



(i.e., iV = 88 — 126) and display a range of ground state 
shapes wide enough so as to be considered a very chal- 
lenging testing ground for a comparison of the (mean 
field) ground state properties predicted with the new 
parametrizations DIN and DIM against the standard 
DIS Gogny functional. Our calculations also included 
the isotopes i60-i64p^ which turn out to be spherical 
and will not be further discussed in the present study. 

The paper is organized as follows. In Sec. |n]we present 
a brief description of the theoretical framework used. The 
results of our study are discussed in Sec. IIIIl There, 
we discuss, in Sec. IIII Al our Gogny-HFB calculations, 
for which axial symmetry is preserved as selfconsistent 
symmetry, used to construct Potential Energy Curves 
(PECs). In a second step we discuss our study of the 
triaxial landscape, providing Potential Energy Surfaces 
(PESs), by constraining on both /3 and 7 quadrupole de- 
formations. In our axial and triaxial HFB calculations 
we considered at the same time, the three parameter sets 
DIS, DIN and DIM. The interaction Gogny-DlS, taken 
as a reference in the present study, is already considered 
as a global force able to describe reasonably well low en- 
ergy experimental data all over the nuclear chart (see, for 
example, [tI. IssI, [Sll. [52l. IsH. [76l - [87{ and references therein) . 
This is also likely to be the situation with the new Gogny 
interactions DIN [H, [zi] and DIM Q but stiU fur- 
ther explorations are required. Therefore, we consider 
all these interactions/functionals in the present study to 
asses to which extent the fine details of our mean field 
predictions for the nuclei i66-204p^ ^^^.g independent of 
the particular version of the Gogny force employed. Af- 
ter discussing the mean field systematics of deformation 
for the considered nuclei, we turn our attention, in Sec. 
IIII Bl to the underlying single-particle properties as func- 
tions of both deformation (axial and triaxial) and mass 
number. This is relevant if one keeps in mind that, from 
a mean field perspective, shape changes arise when the 
deformed single-particle levels are energetically favored 
to a different degree in open shell nuclei (Jahn- Teller ef- 
fect [1^). We also consider the behaviour with neutron 
number of the spherical shell occupancies corresponding 
to the ground state of the different Pt isotopes in order 
to shed some light into the phenomena involved in the 
different deformation regimes. Finally, Sec. IIVI contains 
the concluding remarks and work perspectives. 



II. THEORETICAL FRAMEWORK. 

In order to compute both PECs and PESs we have 
used the (constrained) HFB method together with the 
parametrizations DIS, DIN and DIM of the Gogny in- 
teraction. The solution of the HFB equations, leading to 
the vacuum I$hfb), was based on the so called gradient 
method [i^ll^ to locate the minima. The kinetic energy 
of the center of mass motion has been subtracted from 
the Routhian to be minimized in order to ensure that the 
center of mass is kept at rest. The exchange Coulomb en- 
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FIG. 1: Potential energy curves for the isotopes ifif'-2()4p^ functions of the axial cjuadrupolc moment Q20 calculated with 
the parametrizations DIS (continuous line), DIN (dashed line) and DIM (dotted line) of the Gogny interaction. 
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FIG. 2: (Color online) Q — 7 planes computed with the Gogny-DlS force for the isotopes ^ ~ Pt. The range of Q values 
has been reduced to focus on the interval around the minima. The contour lines extend from the minimum up to 2 MeV 
in steps of 0.25 MeV. Blue (black) contours are the three lowest, green (light gray) ones the next three and red (dark gray) 
contours correspond to the three with higher energies. For each color the full line corresponds to the lower energy, dashed to 
the next contour and dotted to the higher energy. The minimum of the triaxial landscape can be identified by the small ellipse 
surrounding it. 
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ergy was considered in the Slater approximation and we 
neglected the contribution of the Coulomb interaction to 
the pairing field. 

The HFB quasiparticle operators have been expanded 
in an Harmonic Oscillator (HO) basis containing enough 
number of shells (i.e., N= 13 major shells) to grant con- 
vergence for all values of the mass quadrupole operator 
and for all the nuclei studied. Energy contour plots will 
be shown in the ((3,7) plane [H^l (instead of (/3,7)) with 



i20 



HFB 



\2z' 



y^l^HFB) 



\/3 

Q22 = — ($HFB|a:^ - y^l^HFB) 



?20 



^22 



tan 7 



Q2 

Q 



20 



(1) 



(2) 



(3) 



(4) 



Other interesting pieces of information coming from 
the mean field are the single particle energies (SPEs) 
for protons and neutrons. In our calculations, with the 
Gogny interaction, we are solving the HFB equations 
and therefore the only quantities that can be properly 
defined are the quasiparticle energies. However, in or- 
der to have the more usual, Nilsson-like, diagrams we 
have chosen to plot the eigenvalues of the Routhian 

= < + r — A20Q20 — A22(522, with t being the kinetic 
energy operator, F the Hartree-Fock field. The term 
A20Q20 + A22Q22 contains the Lagrange multipliers used 
to enforce the corresponding constraints. We have first 
performed calculations restricted to axially symmetric 
shapes and in a second step triaxiality is included in our 
mean field analysis. In the first case, obviously, the term 
A22Q22 is missing. In addition, the usual mean field con- 
straints on both neutron and proton average numbers 
are taken into account. Parity and time-reversal are self- 
consistent symmetries in our axial calculations whereas 
parity a,nd simplex are the ones imposed in the triaxial 
case 1631. 



III. DISCUSSION OF THE RESULTS. 

In this section, we discuss the results of the present 
study. The systematics of deformation obtained for the 



166-204 



Pt is described in Sec. IIII Al Single par- 



isotopes 

tide properties are considered in Sec. IIII Bl 



A. Mean field systematics of deformation for 

166-204 

The PECs obtained for the isotopes i66-204p^ ^j^j^ 
our constrained Gogny-HFB calculations preserving ax- 
ial symmetry, are shown in Fig. [T] as functions of the 
quadrupole moment Q20- Both prolate (Q20 > 0) and 
oblate (Q20 < 0) sides are displayed. The prolate side 
is equivalent to the triaxial results, to be discussed later 
on, with Q = Q20 and 7 = 0° whereas the oblate side is 
equivalent to Q — IQ20I and 7 — 60°. 

As we can see, the interactions DIN and DIM pro- 
vide PECs which are extremely similar to the ones ob- 
tained with Gogny-DlS. The deformations of the oblate 
and prolate minima are practically independent of the 
force. The axial quadrupole moment Q20 correspond- 
ing to the absolute minima of the PECs increases until 
A w 180 {N w 102) when it reaches the value Q20 ~ 10 b. 
A sudden prolate to oblate shape change occurs around 
A — 188 for all the Gogny interactions considered in the 
present study. Beyond A = 190, absolute oblate minima 
are obtained with quadrupole moments decreasing until 
the spherical shape is reached for ^''"'Pt {N = 126). The 
opposite situation occurs with the secondary minima. On 
the other hand, the depth of both prolate and oblate wells 
(as compared to the spherical maximum) increases with 
increasing N up to A = 182 (roughly mid-shell) and it 
decreases from there on being the decrease more pro- 
nounced for the prolate wells. This, as explained in Ref. 
[63t can be understood first as a consequence of the filling 
of down slopping levels coming from the high j, unique 
parity 113/2 neutron orbital that would explain the in- 
crease of the depth and then, at mid shell, the filling of 
the up slopping levels that would lead to the decrease of 
the height of the wells. 

Slightly lower spherical barrier heights (i.e., the differ- 
ence between the energy of the absolute minimum of the 
PEC and the energy of the spherical configuration) are 
predicted by the new Gogny forces DIN and DIM. Such 
a sensitivity of the spherical barriers with respect to de- 
tails of the effective interactions used has already been 
found in previous studies [1, H, [49, 90]. In our calcula- 
tions, the largest and smallest values of the total pairing 
energies corresponding to the spherical configurations are 
obtained with the sets DIM and DIS, respectively. The 
set DIN provides pairing energies in between. This al- 
ready reflects the different pairing content of the consid- 
ered Gogny functionals but we postpone a discussion on 
this point for later on. 

The absolute values |Ai?p_o| of the energy differences 
between the prolate and oblate minima of the PECs ex- 
hibit two bumps, the first with a maximum at ^ = 178 
(2.19, 1.97 and 2.11 MeV for DIS, DIN and DIM) cor- 
responds to isotopes with a prolate ground state and the 
second one at A = 198 (1.16, 0.95 and 0.99 MeV for DIS, 
DIN and DIM) corresponds to isotopes with an oblate 
ground state. The nucleus ^**Pt, separating the two re- 
gions, has almost degenerate prolate and oblate minima 




FIG. 3: (Color online) The same as Fig. [2] but for the Gogny-DIN force. 




FIG. 4: (Color online) The same as Fig. [2] but for the Gogny-DIM force. 
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FIG. 5: Mean field excitation energies A_E computed with the Gogny interaction DIS (continuous line), DIN (dashed line) 
and DIM (dotted line) are displayed as functions of the deformation parameter 7 for fixed values of the quadrupole moment Q 
corresponding to the lowest minima of the axially symmetric calculations (see Fig. [T]). Results for ^"^Pt are not included due 
to the presence of a spherical ground state. 
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with |A£:p_o| values of 45, 111 and 71 keV for DIS, DIN 
and DIM, respectively. 

The previous axially symmetric results, agree well with 
the ones in Ref. [1] using the parametrization SLy4 f9T| of 
the Skyrme interaction in the particle-hole channel plus a 
zero range and density dependent pairing interaction (92| 
(with strength g = 1000 MeV fm'^ for both protons and 
neutrons) and with previous Skyrme-HF+BCS calcula- 
tions with the parameter set SIII [s^l- They also agree 
well with the results of the axial calculations reported 
in Ref. [6l| using the parametrizations NLl and NL2 
of the relativistic mean field Lagrangian and with ax- 
ial macroscopic-microscopic calculations reported in Ref. 
[s^ j. However, in our axial calculations, prolate and 
oblate minima lie quite close in energy (jAii'p-ol < 2.2 
MeV). Thus, a 7-path connecting them would be possi- 
ble and triaxiality could play a role, as will be discussed 
below, converting some of the axially symmetric minima 
into saddle points. Therefore, in a second step, we have 
also explored the triaxial landscape and construct PESs 
for all the considered nuclei. 

The PESs obtained with the Gogny sets DIS, DIN and 
DIM are shown in Figs. [H [3] and H] in the form of Q — 7 
planes. In order to simplify the presentation, the range 
of Q values plotted is reduced toOb<Q< lib and 
the contour lines are also severely reduced by consider- 
ing contours every 250 keV up to 2 MeV higher than the 
energy of the minimum. We can see, the spherical struc- 
ture in 204pt while 202-i98p^ exhibit oblate (7 = 60°) 
ground states with Q « 3-4.5 b. An island of triaxiality, 
centered around the nucleus ^^®Pt, is clearly visible from 
these figures. The ground state triaxial coordinates ((3,7) 
within such an island of triaxiality evolve from {Q ^5 b, 
7 w 50° - 58°) in ^^^pt to (Q w9 b, 7 w 10°) in ^^^pt. 
The depth of these triaxial ground states, compared with 
the axial ones, is very small (see below). In the case 
of ^^^^^''^Pt, our calculations predict prolate (7 = 0°) 
ground states with deformations Q « 10 b. Again, for 
i76-i72p|- ygj.y shallow triaxial minima are predicted. In 
the case of ^^^Pt, for example, we find {Q ^5 b, 7 « 20°). 
The isotopes i^o-ieep^ display prolate ground states with 
Qr:!3-4.5 b. 

In order to obtain a more quantitative understand- 
ing of the PESs we have plotted in Fig. [Sj the mean 
field energies corresponding to the lowest axial minima 
Q = Q20 of the PECs in Fig. [T] as functions of the de- 
formation parameter 7. We observe that only one of the 
two axial minima remains in most of the cases. The nu- 
clei 202- 198 pt exhibit oblate (7 — 60°) absolute minima 
and the prolate (axial) solutions become saddle points 
with excitation energies AE < 2 MeV. On the other 
hand, i96-i88p-(^ g^j-g j-qxi^q^ 7-soft with triaxial minima al- 
most degenerate with the (axially symmetric) prolate and 
oblate solutions. In the case of -'^^^Pt, for example, we ob- 
tain \AE\triaxiai-saddie < 0.9 McV. Still iusidc the island 
of triaxiality, the oblate configurations in both i86,i84p-(- 
already show the tendency to increase their excitation 
energies. A similar trend for the oblate solutions is ob- 



served within the mass range 182 < A < 174. Oblate and 
prolate configurations for the nuclei i'i'2-i66p|- ^^^.g q^j^g 
close (0.3 < A£: < 1.2 MeV) and softly linked along the 
7 direction. 

A detailed account of the evolution of the ground state 
triaxial coordinates (Q, 7) as functions of the mass num- 
ber A is presented in the left panel of Fig. [6] for the sake 
of completeness. The striking similarity of the ground 
state deformations obtained with the Gogny interactions 
DIS, DIN and DIM becomes evident from this plot. We 
observe the emergence of weakly oblate ground states 
for the isotopes 202-i98p^^ ^^yc other hand, the sud- 
den shape transition observed in the framework of the 
axially symmetric HFB calculations is now replaced by 
a smooth shape change through the island of triaxial- 
ity represented, in our case, by the isotopes i84-i96p^^ 
A prolate deformed regime is predicted for the isotopes 
i78-i82p^^ We find that, the trend of shape changes pre- 
dicted by our calculations in the considered Pt isotopes 
agrees well with the ones obtained in Refs. [sil, HI] and 
the conclusions extracted from the combination of total 
routhian surface calculations plus quasiparticle random 
phase approximation (TRS-f-QRPA) and IBM models in 
Ref. [22] ■ The general trend in our calculations is also 
consistent with results obtained in the framework of the 
Strutinsky approach [H, |35| . Further down in neutron 
number, the TRS results 0, 111] predict a rapid change 
to a triaxial shape also found in our calculations around 
"2pt. Finally, our PESs for "o-ieep^ gxhibit features 
that could be interpreted as the onset of a more pro- 
nounced vibrational character for these nuclei [22j |. 

The rather involved behavior of the neutron and pro- 
ton pairing energies Epp (with opposite sign) versus mass 
number for the Pt isotopes is shown in the middle panel 
of Fig. [S] As expected for a pure TV = 126 shell closure, 
neutron pairing collapses for the nucleus ^"''Pt. We ob- 
serve that proton pairing shows a non constant behavior 
in spite of having constant proton number (Z = 78) that 
comes from selfconsistency effects. On the other hand, 
neutron pairing energies are lower inside the region be- 
tween A = 174 and A — 186 that is precisely where 
the strong prolate deformation develops. The lowering 
of pairing energies is a consequence of the lowering of the 
level density that is needed (Jahn- Teller effect) to induce 
the deformed minima. For other values of A, the neutron 
level density around the Fermi level is higher and, as a 
consequence, pairing correlations are stronger. Concern- 
ing different values of the neutron and proton pairing, 
for the three Gogny functionals the pairing energies fol- 
low the same isotopic trend and the only relevant change 
is in the absolute value that tends to be slightly lower for 
DIS. At this point it is worth to remember that the value 
of the pairing energy shown is related to the amount of 
pairing correlations present in the system but it is by far 
not certain that the correlation is linear, in other words, 
the different values of Epp for different interactions do 
not necessarily imply the same quantitative behavior for 
pairing correlations. 
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FIG. 6: Ground state triaxial coordinates (Q,7) (left panel), neutron and proton pairing energies (middle panel) and Thouless- 
Valatin moments of inertia (right panel) are plotted as functions of the mass number A. Calculations have been performed with 
the Gogny interactions DIS, DIN and DIM. 



Finally, the Thouless-Valatin moments of inertia 
= 3/ (£2+ - Ea+ ) of the first 2+ states in i66-204p^ 
are plotted in the right panel of Fig. [Bias functions of the 
mass number A. They are compared with the experimen- 
tal values extracted from the available systematics for the 
excitation energies of the first 2+ states in i66-204p|. ^ggg^ 
for example, j22|). The energies needed for the computa- 
tion of J^^^ have been obtained using the selfconsistent 
cranking approximation introducing the usual time re- 
versal breaking term —luJx and the subsidiary condition 
(J^) = ^1(1 + 1) [3, [el El- The Thouless-Valatin mo- 
ments of inertia strongly depend on pairing and there- 
fore a comparison of the results obtained with the three 
Gogny functionals considered in this work, which exhibit 
different pairing contents, can also give a hint on the qual- 
ity of their predictions. As can be seen, the results follow 
the same isotopic trend irrespective of the Gogny force 
with the tendency to be the largest for Gogny-DlS and 
the smallest in the case of Gogny-DIM. Nevertheless, the 
differences in the predicted values can still be attached 
to the uncertainties in the effective interactions and we 
observe how the selfconsistent cranking results tend to 
overestimate the experimental values. This defect, well 
known already, is a direct consequence of too low pair- 
ing at the mean field level and its solution would require 
an improved treatment of pairing correlations. There are 
many mechanisms beyond mean field that modify the 
amount of pairing correlations in a given nuclear system, 
but there are two particularly important that tend to in- 
crese correlations. One is the restoration of the particle 
number symmetry broken by the HFB method and the 
other is shape fiuctuations around the HFB minimum. 
The latter is connected with the fact that the HFB mini- 
mum corresponds to a low level density region of the sin- 
gle particle spectrum (see next subsection) and therefore 
its amount of pairing correlations is far lower than the 



one of the neighbouring configurations. Taking into ac- 
count the first mechanism would involve particle number 
projection, whereas the second can only be treated in the 
scope of the GCM with the quadrupole moment as gener- 
ating coordinate or in the Bohr Hamiltonian method [STj . 
Glearly, both methods are out of the scope of the present 
study. On the other hand, we observe certain correlation 
between the evolution with the number of neutrons of the 
quadrupole moments Q and the moments of inertia J^^\ 
as it is apparent by looking at the left and right panels of 
Fig. [HI This correlation is not so evident for the 7 degree 
of freedom. 

The results discussed in this Section indicate, that the 
new interactions/functionals DIN and DIM provide the 
same quality of mean field ground state predictions for 
the considered Pt isotopes as compared with the Gogny- 
DIS force taken as a reference in our calculations. This 
coincidence give us confidence on the robustness of our 
mean field predictions with respect to a change in the 
particular parametrization of the Gogny interaction used. 
The agreement is also rather good with the mean field 
picture obtained in Refs. [1, 111] within the HF-I-BCS 
framework based on the Skyrme parametrization SLy4 
in the particle hole channel plus a zero range and den- 
sity dependent pairing (with strength g — 1000 MeV fm'^ 
for both protons and neutrons). Both Gogny-DlS and 
Skyrme-SLy4 represent well reputed interactions whose 
reasonable predictive power has already been thoroughly 
tested all over the nuclear chart and it is very satisfying to 
observe how the new parametrizations DIN and DIM, in 
spite of the relaxation of some of the original constraints 
in their fitting protocols and more oriented to reproduc- 
ing nuclear masses [H, [tI] , still follow very closely the 
fine details predicted with Gogn y-D lS . Skyrme-SLy4 and 
other theoretical models [1, [2I lli-lll IM HI, for an 
isotopic chain with such a challenging shape evolution. 
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FIG. 7: (Color online) Upper panels: Proton (left panel) and neutron (right panel) SPEs for the nucleus ^^''Pt as functions of 
the axial quadrupole moment Q20- The Fermi levels are also plotted with a thick (red) dotted line. The results correspond to 
the force Gogny DIS. Solid (dashed) lines are used for positive (negative) parity states. With increasing K = 1/2, 3/2, 5/2, . . . 
values color labels are black, red, green, blue, dark-blue, brown, dark-green, etc. Lower panel: the same as above but for 
neutron SPEs of the nuclei ^**Pt (left panel) and ^"^Pt (right panel). The spherical quantum numbers at Q20 ~ are given 
for a number of orbitals close to the Fermi level. 



Let us remark that we are perfectly aware of the fact 
that the (static) mean field picture described above, 
should also be extended to a dynamical treatment of the 
relevant degrees of freedom. This becomes clear from 
the topology of the PESs indicating that, in order to 
access a quantitative comparison of the energy spectra 



and reduced transition probabilities with the considered 
Gogny interactions, the dynamical interplay between the 
zero point motion associated with the restoration of bro- 
ken symmetries (mainly, angular momentum and parti- 
cle number) and fluctuations in the collective parameters 
(/3,7) should be taken into account. For such a cum- 
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FIG. 8: (Color online) Proton SPEs are plotted for the nucleus ^**Pt. In the left panel, the SPEs are plotted as a function of 
the axially symmetric quadrupole moment Q20 from Q20 = up to Q2(i ~ 7.8 b. In the middle part the triaxial proton SPEs 
are plotted as functions of the deformation parameter 7 and for Q — 7.8 b. Finally, on the right panel, the axially symmetric 
SPEs are plotted as a function of Q20 from Q20 ~ —7.8 b up to Q20 ~ 0. The proton Fermi level is depicted as a thick (red) 
dotted line. Some relevant K values are also included. In the axial plots, at Q2Q = (i.e. sphericity) we have the (1^/2 at 
e ~ —9.5 MeV, at an energy of around -7.5 MeV we have the ^3/2 and at e ~ —7.2 MeV we have the /iii/2 orbital. Also at 
Q20 = the Z = 82 shell gap is clearly visible. Solid (dashed) lines are used for positive (negative) parity states. For more 
details, see the main text. 



bersome and computer power deraanding extension the 
Gaussian Overlap Approximation (GOA) appears as a 
first suitable choice [82, 94, 95]. For a very recent and 
excellent pedagogical review, the reader is also referred 
to [9^. Work along these lines is in progress and will be 
reported elsewhere. 



B. Single particle properties 

In this section, we pay attention to single particle prop- 
erties of the considered isotopes. To this end, we first 
show in Fig. [7] SPE plots, as functions of the axial 
quadrupole moment Q2o- The isotopes ^®°Pt, ^*^Pt and 
^°^Pt are taken as illustrative examples. The SPEs cor- 
respond to our Gogny-DlS calculations. For other nuclei 
and/or Gogny interactions the results are quite similar. 

The energy levels in Fig. [7] gather together around 
the spherical configuration Q2q=0 forming the spherical 



shell model orbitals nlj. Due to axial and time- reversal 
symmetries SPE levels, tagged by the K quantum num- 
ber corresponding to the third component of the angular 
momentum in the intrinsic frame, are doubly degenerate. 
Positive and negative parity states are plotted with full 
and dashed lines, respectively. The proton Xz and neu- 
tron Ajv Fermi levels are also shown with a thick (red) 
dotted line. As it is well known, atomic nuclei "avoid" re- 
gions with high single particle level densities (Jahn- Teller 
effect) and therefore the plots of SPEs versus deforma- 
tion help us to identify regions where energy gaps favor 
the appearance of deformed minima 97 1 . 

We also look at the onset of deformation in the con- 
sidered Pt nuclei using the Federman-Pittel (FP) criteria 
(98| . The origin of nuclear deformation is certainly a 
much more complicated phenomenon inv olvin g different 
mechanisms (see, for example, d, [5l|, IH, llOOf and refer- 
ences therein for detailed discussions on this issue) . Even 
though the FP argument may be incomplete, it is cer- 
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FIG. 9: (Color online) The same as Fig. |8]but for neutrons. In the axial plots, at Q20 = (i.e. sphericity) we have the ^19/2 
and /7/2 at e r; —10.5 MeV, at an energy of around -8 MeV we have the 113/2 and at e ~ —7.5 MeV we have the /5/2 and P3/2 
orbitals. Also at Q20 = the N = 126 shell gap is clearly visible. 



tainly playing a role in the onset of nuclear deformation 
and we resort to it in the present mean field study as a 
way to establish a very qualitative and simplified overall 
picture of shape changes in terms of the evolution of the 
und erlying S PEs. For recent studies using similar ideas, 
see [6J, Il0l| . In order to incorporate the terminology of 
spherical orbitals to discuss deformed configurations, we 
assign to a given deformed single particle orbital the label 
of the spherical orbit from which it originates at Q20 = 

In Fig. [71 we can see at zero deformation the spherical 
proton shells 3si/2, 2^3/2 and 2^5/2 below the 

Fermi level Xz and the I/19/2 level above it. The relative 
position of Xz, with respect to the SPEs, is quite stable. 
The main effect of the different neutron numbers appears 
in the scale of Fermi energies Xz ranging from around -9 
MeV in the case of ^°^Pt to values close to zero in the case 
of the neutron deficient isotope ^^^Pt. Therefore, in Fig. 
[7] we have only plotted the proton SPEs for the isotope 
i80p|. rpj_^g closeness to the proton magic number Z = 82 
makes Pt isotopes {Z = 78) to display a tendency to be 
spherical or slightly oblate as it can be seen in Fig. [7]from 
the huge spherical gap of '-^ 6 MeV and the large gap of 



~ 3 MeV on the oblate side centered around Q20 — —5 
b. 

For neutrons the relevant spherical orbitals shown in 
Fig. [7] are I/19/2, 2/7/2, 1^3/2, 2/5/2, 3^3/2, 3pi/2 and 
159/2- Notice also that in the case of ^"''Pt the tendency 
to be spherical is reinforced by its magic neutron number 
N — 126 and the final result is a spherical nucleus. In our 
calculations, the N — 126 spherical shell gap is observed 
to change with mass number from ^ 5.5 MeV in ^*°Pt to 
~ 4.5 MeV in ^^^Pt going through 5 MeV for ^^^Pt. 

Between 202pt and ^^^Pt, Xn crosses a region on the 
oblate side where an energy gap of ~ 3 MeV (i.e. half the 
size of the spherical gap) is found. This occurs at around 
Q2Q = —5 b, which overlaps perfectly with the gap on the 
oblate proton sector already mentioned above. On the 
other hand, as can be seen from Fig. [3 the \iiz/2 gets 
more and more occupied for increasing neutron number 
N and at a certain point (i.e., beyond ^^^Pt) its role is 
transferred to the 2/5/2 and 3p3/2 orbitals. According to 
FP only the neutron 2/5/2 can interact with the proton 
2^3/2 (i.e.. Up = n„, Ip = In — 1) but since the I values 
are in this case low we should not expect a very strong 
interaction among them. This, qualitatively explains the 
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FIG. 10: Spherical occupancies in the proton (left panel) and neutron (right panel) ground state wave functions for the isotopes 
i66-204p^_ For details, see the main text. 



appearance of the, very soft and close to spherical, oblate 
minima displayed in Fig. [1] within the mass range 196 < 
A < 202. A similar mechanism leads to very shallow and 
weakly prolate secondary minima in i96,i98p|-_ 

Between ^^''Pt and ^^''Pt the neutron Fermi level on 
the prolate sector crosses a region with energy gaps of ~ 
3 MeV and between Q20 = 3 b and 10 b. From the neu- 
tron SPE plot of ^^^Pt, we observe how within this mass 
range, the most prominent role is played by the lii3/2 
which, according to the FP criteria, interacts optimally 
with the proton I/111/2 to favor a prolate shape. On the 
oblate side, the interaction between the relevant orbitals, 
leading to the corresponding minima, takes place around 
Q2Q = —5 b. The final net effect of these driving forces 
is the appearance of prolate and oblate minima that be- 
come saddle points inside the island of triaxiality (see 
also the discussion below). 

Below ^*^Pt, our axially symmetric calculations pre- 
dict strong prolate deformations. In this case there is a 
strong gap on the neutron sector centered at Q2a— 10 b, 
as it can be seen from the neutron SPE plot of ^'^'^Pt. The 
most active deformed orbitals are the ones coming from 
the neutron 1ji3/2 and 2/7/2. They interact with the de- 
formed proton orbitals coming from the l/iii/2 and 2^5/2 
to produce strongly prolate deformed shapes. The sec- 
ondary oblate minima predicted for these nuclei, can be 



associated with the proton and neutron energy gaps ob- 
served around Q20— -5 b. Further down in neutron num- 
ber, the neutron Fermi level A at explores regions lower 
in energy resulting in less pronounced prolate minima lo- 
cated around Q20= 4 b in the neutron deficient isotopes 

166,168p^ 

Let us now turn our attention to the origin of triaxiality 
and for this, proton and neutron SPE plots are shown in 
Figs, m and [HI respectively, as functions of the triaxial 
deformation parameter 7 [y, |6^ for the nucleus -'^^^Pt. 
We consider Q =7.8 b which corresponds to a rather 
large region in the Q — 7 plane near the triaxial minimum 
(located at Q —6.2 b, see left panel of FiglHl) where the 
PES is very flat in the two coordinates (see Fig. 
The value of Q =7.8 b also corresponds to the position 
of the axial prolate minimum (see Fig. [T]). These plots 
allow us to identify the K values of the triaxial SPEs at 
the prolate (7 = 0°) and oblate (7 = 60°) limits and 
the change of the K contents observed in most of the 
levels as 7 evolves. Typical examples in this context, 
are the negative parity K=l/2 proton level with SPE 
e « -3 MeV at Q=7.8 b and 7 = 0° which transforms 
into a if = 9/2 level at 7 = 60° and the positive parity 
K = 13/2 neutron level with SPE e w -4.2 MeV at (5=7.8 
b and 7 = 0° which transforms into a, K = 1/2 level at 
7 = 60°. The rather low level density below the proton 
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Fermi level for 7 between 0° and 30° favors the flateness 
of the energy curve as a function of 7 helping thereby the 
development of the triaxial minimum in ^*^Pt around 7 = 
30°. From this plot we also conclude that the isotopes 
with two protons less (Os isotopes) will be more prone 
to triaxiallity as discussed in Ref [63|. Concerning the 
behaviour of the neutron SPEs with 7 depicted on Fig. [S] 
we observe a region of low level density in the interval of 
7 between 0° and 20° that would favor the development 
of triaxiallity. From there on the level density increases 
and therefore, in order to develop a triaxial minimum, 
the system is forced to change its Q deformation to lower 
values (see Fig. [2] for the nucleus under consideration, 
^^^Pt). The removal of two or four neutrons makes the 
level density in the range of 7 between 0° and 20° even 
lower than in the case of -'^^^Pt explaining the 7 deformed 
minima observed in Fig. [5] for the nuclei is^asep^^ 
the other hand, the addition of two or four extra neutrons 
leads to a decrease of the level density in the interval 
between 7 k, 30° and 7 = 60° favoring the appeareance 
of triaxial minima in i90.i92p|. gg^-j- c^j-ves in i94,i96p-j-^ 

as observed in Fig. [5] 

In the spirit of the shell model, it is also interesting to 
compute the spherical occupancies Q,7) of the dif- 
ferent Ij orbitals in the (usually) deformed ground states 
|$hfb(Q,7)) of the Pt isotopes studied in this paper. 
They are given by 

I$hfb(Q,7)>(5) 

n m 

where c^ij^ and Cnijm are the creation and annihilation 
operators of spherical harmonic oscillator orbits charac- 
terized by the quantum numbers n, I, j, m. The sum 
in m is introduced to make the quantity ([5|) invariant 
under changes in orientation (and therefore to represent 
a genuine "spherical" quantity). The sum in the radial 
quantum number n does not allow to pin down which 
specific HO orbital is occupied but, on the other hand, 
allows us to get rid of the uncertainties associated to the 
fact that the radial wave function of the nuclear orbitals 
is close but not exactly the one of the HO. 

The proton spherical occupancies in the ground state 
wave functions of all the Pt isotopes considered are shown 
on the left panel of Fig. [TUl As the number of protons re- 
mains constant along the isotopic chain one could expect 
a flat behavior of all the occupancies. However, we ob- 
serve in that figure how the occupancies of the different 
orbitals can be classified in two different regimes, namely 
the weak deformation regime including i66-i72p^ 
i88-204p|- g^j-^j strong deformation regime including 
the nuclei i74-i86p^ jg noteworthy that in each of the 
regimes the proton occupancies remain rather constant 
irrespective of the 7 deformation and even the specific 
value of the Q deformation parameter. The strong de- 
formation regime differs from the weak deformation one 
in that the d^/2i and /111/2 orbitals loose occupancy 
in favor of the /7/2J fb/2 and hQ/2- This rearrangement 
of the occupancies is mainly due to the smearing out 



of the Fermi surface as a consequence of the increasing 
proton pairing correlations (see Fig. ^ and to a lesser 
extent to the quadrupole interaction among orbits that 
can transfer particles from one orbit to another. Given 
the shift of occupancies from the ^5/2, ^3/2 and /111/2 to 
the /y /2 , /s /2 and /19 /2 orbitals we can interpret the well 
deformed ground state of i74-i86p-j. ^ multiparticle- 
multihole excitation out of a reference spherical ground 
state. The number of particles exchanged in this kind of 
spherical shell model language is in between two to four 
according to the results shown in the left panel of Fig. 

m 

In the case of neutrons, shown on the right panel of 
Fig. 1101 as neutron number increases we are occupying 
orbitals belonging to the N = 5 negative parity major 
shell and the positive parity intruder ii3/2- As a conse- 
quence of deformation, the spherical orbitals are mixed 
up and therefore placing two particles in a given deformed 
orbital by no means implies placing two particles in an 
spherical orbit. The two particles will be distributed 
among all the components of the deformed orbital when 
expressed in the spherical basis. As a consequence, the 
behavior of the spherical occupancies with neutron num- 
ber is more or less linear in the whole interval and for 
all the orbitals involved. The only noticeable deviation 
from this trend takes place when entering the strong de- 
formation regime where the /7/2, ^19/2 and /111/2 orbitals 
loose particles in favor of the high j orbitals 113/2, *ii/2j 
ji3/2 and J15/2. Also the 59/2 orbital gets more parti- 
cles (through the coupling to the = 6 orbital). This 
change in occupancies can be mostly attributed to the 
quenching of neutron pairing correlations in the ground 
state of the strongly deformed isotopes i74-i86p^^ When 
the weak deformation regime is entered at ^ = 188 a 
much smoother behavior of the spherical occupancies is 
recovered. This, together with the smooth behavior of 
proton occupancies is quite unexpected as in the weak 
deformation regime the different isotopes have a variety 
of ground state deformations ranging from prolate to tri- 
axial to oblate. The conclusion is that the occupancies 
are more sensitive to the magnitude of the deformation 
Q than to the 7 degree of freedom. 

The proton levels closest to the Fermi level are the 
■Si/2j ^ii/2j £^3/2 and c?5/2 and therefore are the ones ex- 
pected to strongly interact with the neutron spin orbit 
partner according to the FP mechanism. The more ef- 
fective orbitals are those with high j and therefore the 
proton /111/2 is expected to strongly interact with neu- 
trons hg^2 which is empty at the beginning of the chain 
and gets steadily occupied as more neutrons are added. 
Also a strong interaction is expected with the neutron 
*i3/2 orbital that shares with the /ig/2 one the occupa- 
tion pattern as a function of N. It is also noteworthy to 
point out how proton hg/2, that should be empty accord- 
ing to the spherical shell model, gets some occupancy for 
the well deformed nuclei i''4-i86pt consequence of 
the smearing out of the Fermi surface, that allows those 
orbitals to interact via FP with the neutrons ft.11/2 and 
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(^9/2. It could be argued that the proton /111/2 orbital 
is loosing two particles in the strong deformation regime 
and this could imply a strong impact in its interaction 
with its neutron spin-orbit partner. This would be true 
if the ^11/2 orbital were occupied by a small amount of 
neutrons (of the order of two) as this would imply com- 
pletely emptying out the orbtial, but it has to be bear in 
mind that there are ten neutrons at the begining of the 
region under consideration. 



IV. CONCLUSIONS 

In this paper we have studied the evolution of the 
ground state nuclear shapes in a series of Pt isotopes 
ranging from N = 88 {A = 166) up to TV = 126 {A = 
204), covering practically one complete major shell. The 
study has been performed within the selfconsistent HFB 
approach based on the DIS fs^ and the recent DIN [t^I 
and DIM 75] parametrizations of the Gogny interaction 
and we have included, in addition to the axially symmet- 
ric limit, the triaxial degrees of freedom /3 and 7. 

From the analysis of the axially symmetric limit, we 
conclude that a sudden prolate to oblate shape change 
occurs around N = 110 (A = 188) for all the Gogny 
parametrizations considered. This result is also in agree- 
ment with those obtained either with Skyrme forces or 
from relativistic mean field calculations ^ HO, [HH [63l |. 
On the other hand, when triaxiality is taken into ac- 
count, the picture that finally emerges is that of smooth 
transitions between the different shape regimes. We find 
that the absolute minimum of the PESs for the Pt iso- 
topes evolves from prolate shapes with increasing values 
of their quadrupole moments Q in the lighter isotopes 
A = 166-182 (with the exception of A = 172-176, which 
exhibit a tendency to triaxiality), to triaxial 7-soft in the 
intermediate isotopes with A = 184 — 196, and to oblate 
shapes in the most neutron-rich isotopes A = 198 — 202. 
Finally, the isotope ^'^''Pt becomes spherical. By ana- 
lyzing PECs and the Q — ^ landscapes, we observe that 
the (axial) prolate and oblate minima, well separated by 
high energy barriers in the /3 degree of freedom, are softly 
linked along the 7 direction. Indeed, most of the sec- 
ondary axial minima become saddle points when the 7 
degree of freedom is included in the analysis. Pairing 
energies and Thouless-Valatin moments of inertia have 
also been analyzed as functions of the number of neu- 
trons, finding some correlation with the evolution of the 
quadrupole deformation Q. Such a correlation, on the 
other hand, is not so evident in the case of the 7 angle. 

We consider the similarity between the mean field 
ground state properties predicted for the considered Pt 
chain with the most recent versions of the Gogny interac- 
tion (i.e., DIS, DIN and DIM) employed in the present 
study as very positive. On one hand, the results give us 
confidence concerning the robustness of the predictions 
against the details of the particular effective Gogny in- 
teraction used. The robustness is reinforced when the 



trend of our calculations is compared with the ones ob- 
tained with Skyrme forces [1, [6J| and other theoretical 
approaches [3, IH, [s^, ll^, which are again very similar. 
On the other hand, our results also point to the fact that 
the new incarnations DIN and DIM of the Gogny inter- 
action, based on fitting protocols more in the direction 
of astrophysical applications [73, 'ts'] , essentially keep the 
predictive power of the Gogny-DlS force 52] already con- 
sidered as a (standard) global force. 

We have analyzed the proton and neutron SPEs as 
functions of both axial and triaxial deformation parame- 
ters in some illustrative examples. The analysis has been 
done in terms of the density of levels around the Fermi 
surface (Jahn- Teller effect [8^|) and the Federman-Pittel 
mechanism [98] . Our discussion has also been illustrated 
with the calculation of the spherical occupancies in the 
ground state wave functions of the considered Pt nuclei. 
As a result, we obtain a qualitative understanding of the 
emergence of deformed configurations with two main in- 
gredients which are, the energy gaps that appear at dif- 
ferent deformations in the SPEs of neutrons when the 
Fermi level Xn crosses different regions, and the special 
role of the overlap between the proton I/111/2 and the 
neutron 1^13/2 orbitals. 

The study of the low- lying excitation spectra and tran- 
sition rates in conjunction with shape transitions, would 
require to extend the present mean field approach to take 
into account correlations related to restoration of broken 
symmetries and fluctuations of collective variables. The 
restoration of broken symmetries would imply, among 
others, triaxial projections to restore rotational symme- 
try and this is a very difficult and delicated issue with 
realistic forces [102-105j. These difficulties lead to con- 
sider instead of the exact projection some kind of approx- 
imation to it that could eventually end up in a kind of 
collective Bohr Hamiltonian [82, 94 96] with deformation 
dependent parameters. Upon completion of this work, a 
preprint has appeared [871] dealing with the calculation of 
2+ excitation energies in the framework of the 5-D Bohr 
hamiltonian method with parameters extracted from a 
microscopic mean field calculation with the Gogny DIS 
force. As a consequence of the number of nuclei con- 
sidered (around two thousand) in that calculation, their 
analysis of the mean field results is not as exahustive as 
ours. 
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